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Abstract 

Direct photon production data at fixed target and collider experiments from pp and 
pp collisions are analysed in NLO QCD. The results are grouped into three sets having 
a/s = 20 — 30 GeV, 630 GeV and 1.8 TeV. It is seen that at the lowest energies considered, 
the theory systematically underpredicts the results for all the values of the transverse mo- 
mentum. At the intermediate energy of the UA2 experiment, the data are more accurately 
described at high pt-, and at the highest collider energies no significant underprediction 
is seen even at low p^, when using recent structure functions as measured at HERA. 
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Experimental and theoretical studies of direct or prompt photon production in hadron- 
hadron collisions at large transverse momenta have been pursued extensively in the lit- 
erature. This interest stems from the expectation that direct photons may provide a 
clean test of perturbative QCD, and also give valuable information about the parton 
distribution functions of the hadrons. The production of direct photons in lowest order 
proceeds through the annihilation {qq —>■ (77) and the Compton {qg — > ^7) processes. 
The Compton processes have the unique distinction of being sensitive to the gluon struc- 
ture function. At the same time a fairly precise next to leading order QCD treatment is 
available for the analysis |T|]. 

This offers us an opportunity to perform a global test of the QCD theories for the direct 
photon production. Huston et al. have recently performed this task using CTEQ2M 
structure functions, and reported that in general the theory underpredicts the data. They 
find a general trend of a larger deviation at lower pt, and a smaller deviation at larger 
Pt at all the energies, which they claim should require an additional nonperturbative pt 
broadening of the initial state parton distribution. 

In this note, we will demonstrate how this apparent gap between theory and data is 
strongly reduced when using more recent structure functions, and how NLO QCD nicely 
describes the data in the region where this is to be expected. 

This analysis relies on the NLO QCD calculation of direct photon production by 
P.Aurenche et al. This treatment was recently applied to the data for direct photon 
production from the Fermilab collider, see for instance [Q, to the global analysis mentioned 
above 0, as well as to an overview of existing data in order to form a basis for the 
prediction of direct photon production at future (LHC,RHIC) experiments 0. 

In a spirit similar to that of ref. we use the NLO QCD theory but with a 

different parametrization of the structure functions, we now compute the cross section 
of direct photon production in pp and pp reactions and compare our results with the 
existing data. The available experimental results for direct photon production in pp and 
pp collisions can be broadly categorized into three groups, depending on the centre of mass 
energy; ^/s = 20-30 GeV, 630 GeV, and 1.8 TeV. This increase of energy corresponds to a 
decrease of the momentum fraction x which is probed by the photon production process 
inside the hadron. The elementary parton-parton collision occurs at s = X1X2S, where 
^/s is the cm. energy of the hadron-hadron reaction, and Xi and X2 are the momentum 
fractions of the parton from the beam and target hadron, respectively. For a photon 
emitted at 6' = 90° in the c.m.s., corresponding to a rapidity oiy = 0, the expressions are 
symmetric in xi and X2, which we then simply denote by x (kinematically, the process 
actually converts a quark into a photon momentum). All of the data we consider here 
are taken around midrapidity. Calculating the production process involves an integration 
over Xi, X2 within the limits Xmin ^ x <1. For y = 0, and using the normalized transverse 
momentum variable xt = 2pT / a/s, we have 



Therefore, the data for the first group should be sensitive to parton structure functions 
down to Xmin=0.16, whereas for the second group this sensitivity should be down to 
a^min=0.02. For the highest energy considered, the structure function is explored down to 



X- 





Pt 




x=0.0067. 
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In our calculation, we use the structure functions MRSD—' from Martins, Stirling 
and Roberts [^], which provides a good fit to recent HERA data. We first illustrate the 
difference to the CTEQ2M structure function, which was used previously in this context, 
by comparing both parametrizations in the x and range which is relevant to the data 
we are considering. Fig. 1 shows the parton distribution functions provided by CTEQ2M 
and MRSD-' for values of = (lo GeV)^ and = (100 GeV)^, which is about the 
range covered by the UA2 and Tevatron photon data. For the gluon distribution, plotted 
in fig. la, one observes a strong effect of the evolution on both structure functions. 
However, the parametrizations are also different. With decreasing x, which corresponds to 
going to smaller values of px, more gluons are provided by MRSD—' than by CTEQ2M. 
Going to higher ^/s, corresponding to probing of smaller x, similarly opens a gap between 
the gluon luminosities provided by the two structure functions. Since the low pt region is 
dominated by the Compton processes involving gluonic fusion, the calculation which uses 
the CTEQ2M structure function therefore results in a lower cross section than the one 
using MRSD—' . For a comparison, we also show the gluon density provided by the latest 
(1995) parametrization MRSG. It has the same trend of providing more gluon luminosity, 
even somewhat more pronounced than MRSD—'. For a better comparison to existing 
calculations, we are going to use MRSD—' in the following calculation. In fig. lb we plot 
the valence {u only) and the sea quark distributions. Here, less difference is seen when 
comparing the two structure functions CTEQ2M and MRSD—' which is of the order 
of the effect of the evolution. Correspondingly, the result for the photoproduction 
cross section for both distributions shows only little difference. This applies to the high 
Pt region where qq annihilation is the dominating process. 

We now turn to the results for direct photoproduction in pp and pp reactions. For the 
three groups of energies, we show the comparison of the NLO QCD calculation and the 
data in the form of the fractional difference {data — theory) /theory using the invariant 
cross section 

da 1 da 

rp /oN 

d^p 271 Pt dptdy 

The results are shown for the different energy regions with increasing energy in fig.'s la-lc. 
The low energy group, shown in fig. la, contains the data of NA24 [Q, WA70 UA6 



and E706 |T^] as compared to our result. Here and in the following, we extrapolated 
data taken in a finite y interval to y = by using the procedure of Aurenche and Whalley 
| |11| . This calculation as well as the following ones are carried out at a fixed scale of 

= Pt/4: as is that of Huston et al. 0]. As discussed for instance in 0, the dependence 
on the scale is weak in this domain and mainly affects the overall normalization, rather 
than the slope. We also note at this point that a possible overall uncertainty in the 
normalization of the data will also not affect the slope, as is discussed in p[. 

What we find in the comparison with the low energy data is a general trend of an 
underestimate of the data by the NLO calculation, however no significant pt dependence 
of the discrepancy. 



At intermediate energy, only the data of UA2 [|I2[ exist and are shown in the com- 
parison to theory in fig. 2b. As is seen in the plot, the data are already better described 
by theory, and the trend of a better agreement at high pj- might be read off. However, 
the absolute deviation is reduced as compared to the calculation based on the CTEQ2M 
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distributions (apparently, the first point in fig. 2b is not contained in fig. 4 of or is 
outside of the scale). 

In the last figure, 2c, we show our results compared to the collider data from DO 
|T3| and CDF |jT^. Here, the agreement is satisfactory, and no significant trend of a 



discrepancy increasing with decreasing pt can be read off. The apparent discrepancy 
between our result and the corresponding one of Huston et al. has its origin in the lower 
gluon luminosity provided by the CTEQ2M parametrization in the x region relevant 
for the low px data, as we discussed before. It is also felt that the introduction of an 
additional broadening of the initial state parton distribution by some value < A/c^ >~ 1 
GeV of presumably nonperturbative origin would not have provided a resolution of this 
discrepancy. For the UA2 as well as for the high energy data, < Akx ><^ Pt at which the 
data are taken and thus the effect of the additional < Akx > broadening would be hardly 
visible. For the data in the low energy region, on the other hand, the underestimate of 
the data is present for all values of px- 

We conclude that by using the recent structure function MRSD—' (and even more 
so by the latest parametrization MRSG), we find that deviations are in general much 
less as compared to the results of Ref. p[, and follow the logical sequence of being large 
at lower energies and smaller px and becoming insignificant at the highest energy over 
the entire range of px, for which the data are available. This conclusion is corroborated 
by the differences in the structure functions used at lower x values. The behavior we 
find is precisely what is expected from any NLO calculation in general, here (9(aa^). At 
small ^/s, the NLO QCD calculation necessarily underpredicts the data, then starts to 
give a quantitative prediction at sufficiently high y/s first at high values of pt, and finally 
provides a satisfactory quantitative description of direct photon production over the entire 
range of pt range at high -^i. 
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FIGURE CAPTIONS 



1.) Comparison of the two structure functions used in the QCD analysis, MRS D-' and 
CTEQ2M, in the x and range relevant to the collider experiment, Q'^ ~ (Pt)^- 
Also shown is the most recent parametrization MRS G. Fig. la: Density distribution 
of gluons, fig. lb: Valence and sea quark distribution. 

2a.) Fractional difference between the NLO QCD calculation and data taken at low 
energies, a/s < 30.6 GeV. 

2b.) Fractional difference between the NLO QCD calculation and data taken at the 
intermediate energy of y/s = 630 GeV by UA2. 

2c.) Fractional difference between the NLO QCD calculation and data taken at the 
Tevatron at ^/s = 1.8 TeV. 



6 



X G(X) 




X S(X), X Uv(X) 




(Data -NLO QCD)/NLO QCD 



[NO 



O 

CO 



O 



o 

P o 



o 

cn 



o 

05 



"1 I I I I I I I r 



^ □ O X 

H a ^ 

^ > > > 

o 05 ^ ro 

05 O 4^ 



I ^ 1 



I <> 1 



J \ \ \ \ \ \ \ L 



"1 I I r 



I 1 



I 4 1 

I 4 1 

I I ^>< I I 



"1 I I r 

a 

1^ 

II 



-3- 



-3- 



\ ^ 



-3- 



I I 
I X 1 



I 1 



-B- 



-3- 



J \ \ \ \ \ \ \ L 



(Data -NLO QCD)/NLO QCD 



J\3 



[NO 



O 

cn 



O o 



O 



o 

h- »• 
cn 



o 



"1 I I I I I I I r 



I — )^ — I 



I — « — I 



I — n — I 
« — I 



I — ^ 



I — « — I 



I ^ 



J \ \ \ \ \ \ \ L 



"1 I I I I I I I r 



I « 1 



I K I 

I K I 
» I 



I « 1 



J \ \ \ \ \ \ \ L 



(Data -NLO QCD)/NLO QCD 



O 
cn 



O 



o 

^ o 



cn 



o 

cn 

"1 I I I I I I I r 



o 
b 



f\3 



cn 



I 



^ 

I 




I — 



I X I 



I ^ 



I ^ 



J ^ \ \ \ \ \ \ L 



o 

CJl 

"1 I I I I I I I r 



n 



a 

o 



J \ \ \ \ \ \ \ L 



